Combining ability analysis of yield and biomass allocation related traits in newly developed wheat populations

Increasing biomass allocation to the root system may increase soil-organic carbon stocks and confer drought adaptation in water-limited environments. Understanding the genetic bases and inheritance of biomass allocation is fundamental for drought tolerance breeding and soil health. The objective of this study was to determine the general and specific combining ability, maternal effects and the mode of gene action controlling the major yield and biomass allocation related traits in wheat to identify good combiners for breeding and enhanced carbon sequestration. Ten selected wheat genotypes were crossed in a full diallel mating design, and 90 F2 families were generated and evaluated in the field and greenhouse under drought-stressed and non-stressed conditions. Significant differences were recorded among the tested families revealing substantial variation for plant height (PH), kernels per spike (KPS), root biomass (RB), shoot biomass (SB), total plant biomass (PB) and grain yield (GY). Additive gene effects conditioned PH, SB, PB and GY under drought, suggesting the polygenic inheritance for drought tolerance. Strong maternal and reciprocal genetic effects were recorded for RB across the testing sites under drought-stressed conditions. Line BW162 had high yield and biomass production and can be used to transfer favourable genes to its progeny. The parental line LM75 maintained the general combining ability (GCA) effects in a positive and desirable direction for SB, PB and GY. Early generation selection using PH, SB, PB and GY will improve drought tolerance by exploiting additive gene action under drought conditions. Higher RB production may be maintained by a positive selection of male and female parents to capture the significant maternal and reciprocal effects found in this study.

www.nature.com/scientificreports/ variation in wheat biomass allocation resulting in different genotype responses to drought stress 14,15 . Hence, promising genotypes can be identified and used in crossing programs to develop breeding populations to select transgressive segregants with high root biomass 15,16 . Increasing the root biomass of cultivated crops will further contribute carbon into soil via root biomass deposits and rhizodeposition that will lead to a net increase in soil carbon stocks for climate change mitigation 17,18 . Nevertheless, focusing on increasing root biomass alone without maximizing yield-related traits may lead to a loss of wheat productivity. In wheat breeding, plant height, kernels per spike and harvest index are key yield-related traits that have been targeted for drought tolerance breeding 19,20 . Multiple traits selection for high biomass production and grain yield-related attributes may be required to increase genetic gains in wheat breeding programs 21 .
Breeding gains can be achieved by understanding the genetic basis and inheritance of yield components and creating desirable progenies through combining ability (CA) analysis 22 . The goals of CA analysis are to identify genetically superior lines with high breeding values and to identify desirable cross combinations to improve average performances for cultivar development 23,24 . Combining ability analysis can be done at F 2 generation without substantial loss of information on the breeding values of parental lines 25 . Several studies have conducted genetic analysis on F 2 populations in self-pollinating crops such as wheat 26,27 , soybean 28,29 , common bean 30 and groundnut 31 .
Estimates of combining ability effects and the extent of variance components reveal the magnitude of both additive and non-additive gene action 32,33 . In hybrid breeding, additive, dominance and epistatic interactions of non-allelic genes influence maximum heterosis 34 . On the other hand, additive gene effects are more important in line breeding with minimal contribution from non-additive gene effects, which are lost during segregation in early generations 35 . As a result, maternal effects (inheritance of cytoplasmic genes from mitochondria and chloroplasts and their interaction with nuclear genes), are often ignored in explaining variation among genotypes 36,37 . In pure line cultivar development, maternal effects can be exploited to identify male and female lines in crosses to maximise genetic gain for traits with significant reciprocal effects 38 .
Diallel analysis can aid in partitioning the general combing ability (GCA) and specific combining ability (SCA) effects. It also reveals the magnitude of maternal effect that will be useful in the breeding of wheat. Maternal effects contribute to early seedling development and biomass production of plants 39 . Understanding the magnitude of variation attributed to cytoplasmic DNA would greatly enhance selection efficiency, including for root biomass. Several studies have reported significant maternal effects in crops for agronomic traits 40 pointed out a maternal effect influencing root traits in F 1 reciprocal crosses in common bean (Phaseolus vulgaris L.). Additionaly 22 , reported that salt stress tolerance levels in reciprocal crosses of wheat were related to the maternal plants of the respective progenies. In a study to improve pre-harvest sprouting in barley 41 , reported that seed dormancy was maternally inherited. Conversely, non-significant maternal effects have been reported for wheat grain yield, though small numbers of crosses reflected significant reciprocal effects 42 . Based on the reports mentioned above, there may be significant maternal and reciprocal effects for root biomass allocation and other traits of interest in wheat.
Biomass allocation is an important attribute in developing drought-tolerant crop ideotypes that contribute to yield gains and carbon sequestration. To our knowledge, there is not enough information on genetic analysis of biomass allocation for trait integration and to guide selection and breeding in wheat. To examine the pattern of genetic inheritance of biomass allocation, 10 selected bread wheat lines with contrasting root biomass and drought tolerance were crossed using a full-diallel mating design, and the developed crosses were advanced to the F 2 generation. Therefore, the objective of this study was to determine the general and specific combining ability, maternal effects and the mode of gene action controlling the major yield-related traits and biomass allocation in wheat to identify good combiners for breeding and enhanced carbon sequestration. Information from this study will help determine the suitable parental selection criteria for the efficient breeding of drought-tolerant wheat cultivars.

Results
Analysis of variance. The combined analysis of variance with means squares and significant tests for biomass traits and yield components for parental lines and direct and reciprocal crosses at the F 2 generation is presented in Table 1. The effects of family and family × site interaction were significantly different for PH, SB, RB, PB and GY but were non-significant for HI. The interaction of family and water regime were only significant for GY. Non-significant differences were observed for family × water regime × site interaction.
Mean performance of parental lines and F 2 families. Drought stress had a substantial impact on all measured traits, with mean PH being reduced by 9.38 cm and 5.94 cm in the field and greenhouse, respectively (Tables 2 and 3). The percentage drop in KPS due to drought stress was 19.18% in the field which was double that of 8.8% in the greenhouse. Drought-stress positively impacted root growth, which increased by 78% in the greenhouse but had a reduced effect on RB in the field. Biomass traits were severely reduced by drought with SB experiencing losses of 41.09 g m −2 and 100.82 g m −2 while PB had losses of 181.40 g m −2 and 327.39 g m −2 in the field and greenhouse, respectively. Grain yield reductions due to drought stress were up to 120.67 g m −2 in the field and 410.31 g m −2 in the greenhouse. This represented a yield loss of 40.27% and 80.21% in the field and greenhouse, respectively. Plant height and KPS had higher mean values in the field than in the greenhouse under both water regimes. On the other hand, higher mean values were observed in the greenhouse than in the field for RB and PB under drought-stressed and non-stressed conditions. Some parents and F 2 families were more drought tolerant than others. For instance, parental lines LM75 and BW162 had the highest grain yield and biomass production than other parents under drought-stressed conditions (Supplementary Table 1). Kernels per spike were high in lines LM26 and LM75 under the same conditions. Combining ability analysis across sites. Parental lines had significant (P < 0.05) GCA effects for all recorded traits (Table 5). Significant SCA effects among F 2 families were observed for PH, KPS, SB, RB and PB. Reciprocal effects were significant for KPS and biomass traits such as SB, RB and PB with reciprocal crosses showing significant maternal effects only for KPS. The GCA × site interaction effect were significant for parents for PH, KPS, SB and PB. Similarly, SCA × site effects for F 2 families were significant for the same traits in addition to RB and GY. Reciprocal effects had significant interaction with sites for all traits except KPS. Maternal effects and site interaction influenced SB, RB and GY.
General combining ability effects. The general combining ability of parental genotypes are recorded in Table 6 for drought-stressed and non-stressed conditions. Under drought stress, parent BW141 had negative GCA effects for PH at both sites and positive GCA effects on GY in the greenhouse condition. Parental line LM26 showed significant and positive GCA effects for KPS with negative GCA's for lines BW140 and LM47. No parents showed significant GCA effects for RB under drought. The GCA effects for SB, PB and GY were significant in a desirable direction for LM75 in the field. Parental lines LM70 had positive GCA effects on PH, SB and PB in the greenhouse. Under non-stressed conditions, BW140 maintained negative GCA effects for PH and PB while LM47 showed significant positive effects for PH in the field. The GCA effects for KPS were strong and positive for LM48 at both sites. Negative GCA effects for RB were observed in the field for lines LM48 and LM70. Genotype LM75 showed positive GCA effects for RB in the field and PH, SB, PB and GY in the greenhouse.
Specific combining ability and reciprocal effects. The SCA effects of direct crosses are shown in Table S3. Similar to the GCA effects of parental lines, no F 2 families maintained positive SCA effects for all the recorded traits across the sites. Positive SCA effects were observed for BW141 × BW152 for PH under drought Table 1. Mean squares and significant tests from a combined analysis of variance of ten bread wheat parental lines, 45 direct crosses and 45 reciprocal crosses in the F 2 generation for yield components and biomass traits evaluated at two sites under drought-stressed and non-stressed conditions. *Significant at P < 0.05; **P < 0.01, ***P < 0.001, number of sampled plants = 5, SOV , source of variation; d.f, degrees of freedom; PH, plant height (cm); KPS , kernels per spike; SB, shoot biomass (g m −2 ); RB, root biomass (g m −2 ); PB, total plant biomass (g m −2 ); GY, grain yield (g m −2 ) and HI, harvest index (%). www.nature.com/scientificreports/ were observed for BW162 × LM48 under non-stressed conditions. Strong reciprocal effects were recorded in the F 2 families LM47 × BW152, LM26 × BW140 and LM75 × LM47 under drought-stressed conditions (Table S4). Table 3. Mean values for yield components and biomass traits of the top two yielding bread wheat parental lines and their F 2 generation of direct and reciprocal crosses evaluated at two sites under non-stressed conditions. Number of sampled plants = 5, PH, plant height (cm); KPS, kernels per spike; SB , shoot biomass (g m −2 ); RB, root biomass (g m −2 ); PB, total plant biomass (g m −2 ) and GY, grain yield (g m −2 ); HI, harvest index (%); GH, greenhouse; LSD, least significant difference; SEM, standard error of mean; CV, coefficient of variance. www.nature.com/scientificreports/   www.nature.com/scientificreports/ Table 5. Summary mean squares and significant tests of combining ability and maternal effects for yield components and biomass traits for a full diallel cross at the F 2 generation evaluated at two sites under drought-stressed and non-stressed conditions. *Significant at P < 0.05; **P < 0.01, ***P < 0.001, number of sampled plants = 5, SOV, source of variation; df , degrees of freedom; PH, plant height (cm); KPS, kernels per spike; SB , shoot biomass (g m −2 ); RB, root biomass (g m −2 ); PB, total plant biomass (g m −2 ); GY, grain yield (g m −2 ) and HI, harvest index (%); GCA, general combining ability; SCA, specific combining ability effects; REC, reciprocal effects; MAT, maternal effects; NMAT, non-maternal effects.  Table 6. Estimates of general combining ability effects of 10 bread wheat parental lines, involving progeny at the F 2 generation for yield components and biomass traits evaluated under drought-stressed and non-stressed conditions in the field and greenhouse. *Significant at P < 0.05; **P < 0.01, ***P < 0.001, number of sampled plants = 5, PH , plant height (cm); KPS, kernels per spike; SB, shoot biomass (g m −2 ); RB, root biomass (g m −2 ); PB , total plant biomass (g m −2 ), GY, grain yield (g m −2 ); HI, harvest index (%); GH, greenhouse.

Discussion
Significant differences were observed among the tested F 2 families (Table 1), revealing substantial variation for PH, RB, SB, PB and GY. This offers an opportunity for an effective selection of biomass-related traits to improve drought tolerance and enhance carbon sequestration of cultivars for sustainable wheat production. The significant interaction of families by site (Table 1) showed that genotype ranking changed in the different sites indicating strong genotype by environment interaction affecting all the measured traits. Among all the assessed traits, GY was highly influenced by the interaction effect of families and water regime (Table 1). This suggests that GY is more sensitive to moisture fluctuations than the other traits, which may affect selection response. Therefore, to achieve yield stability, genotypes should be evaluated across multiple growing environments with different moisture availability to identify stable and drought-tolerant genotypes for water-limited environments 48 . All other traits, including PH, KPS, SB, PB and GY were severely and negatively affected by drought stress (Tables 2 and 3). However, individual genotypes responded to drought stress differently, with some parents and families showing high levels of drought tolerance compared to others. The different phenotypic performances observed allow for the targeted selection of better-performing families for genetic advancement, while parental lines can be selected to develop new breeding populations for either dryland or irrigated wheat production 49 . Some families, including BW141 × LM26, LM47 × BW152 and LM26 × BW140 scored higher than all the parental genotypes for RB, PB and GY, respectively, indicating that parental genotypes were able to transmit favourable genes to their progeny under contrasting levels of soil moisture availability.
Evaluation of the test genotypes was done at contrasting growing conditions, namely greenhouse and field environments. Data on the assessment of families in the greenhouse using pot experiments were included to capture and assess total root biomass in the pots, which cannot be completely recovered in field experiments. In addition, root evaluation in pot experiments and growth chambers is easier, and more accurate as there is no mixture of roots from adjacent plants or different genotypes when collecting root samples 50 . In this study, plants allocated more biomass belowground to promote denser and deeper root growth for efficient water and nutrient uptake under drought-stressed conditions 51 . However, this increase in RB production was high in the greenhouse and small in the field, indicating greater accuracy of greenhouse trials in studies involving root phenotyping. Different drought-induced root growth responses have been observed between greenhouse and field conditions under different stress severity 11 . Many confounding environmental effects in field trials compared to greenhouse trials may contribute to the low accuracy of root phenotyping in the field 52 . Pot experiments represent suitable environments for controlled root evaluations but extrapolations based on such data sets would not accurately represent field trials.
The significance of maternal effects for SB under non-stressed conditions (Table 4) showed that a portion of the non-additive gene action in the population could be captured by choice of parent used as a male or female in the mating design. The GCA/SCA ratio was slightly higher under drought-stressed conditions than under nonstressed conditions showing that water availability was able to alter the proportion of additive and non-additive gene action affecting the traits of interest. Similar results were observed by 48 for soybean. The GCA/SCA ratio was > 0.5 for PH, SB, PB and GY under drought-stressed conditions. These traits can be improved by exploiting additive gene action at early generations, thereby increasing selection efficiency and reducing the breeding cycles 33 . Several studies have reported quantitative trait loci (QTL) with additive gene effects for GY and related traits, allowing for the selection of transgressive progenies by crossing superior parents 33,53,54 . Conversely, under non-stressed conditions, there was a greater influence of non-additive gene action, which opens an opportunity to exploit dominance and epistasis for hybrid breeding in irrigated environments. Heritability estimates were low for all measured traits, indicating that the phenotype was a poor measure of the genetic merit of the evaluated lines and families, which reduces the effectiveness of selection. Similarly, low heritability values have been reported in other studies 55,56 .
The presence of interaction effects of sites with combining ability and maternal effects for the majority of measured traits (Table 5) highlights the confounding effects of the environment on gene expression and the importance of choosing the correct selection environment to assess the genotypic value of genotypes. The environments showed that RB had significant reciprocal effects under drought-stressed conditions across sites. Based on those observations, choosing a male or female parent is important to ensure the inheritance of maternal genes for improved root biomass under drought-stress. Strong reciprocal effects influencing root traits have been reported in wheat for different abiotic stresses such as salt stress 22 and cold stress 57 . Also, in a study evaluating interspecific hybrids of sunflower 58 , reported reciprocal effects on root traits that extended beyond the seedling stage and were expressed in mature individuals. Under non-stressed conditions, the maternal effects were observed for PH, SB, RB, PB and GY but were inconsistent across different sites limiting the usefulness of these effects for breeding purposes.
The GCA effects of parents change from positive to negative for all recorded traits across the test conditions (Table 6). This suggests that the environment can influence gene expression in individual genotypes. Parental line LM75 maintained GCA effects in a positive direction for SB, RB, PB and GY. It was the only genotype with a positive effect observed for root biomass under field conditions. This genotype can be utilized to transmit additive quantitative trait loci for the improvement of biomass production and overall yield potential. Different genotypes were responsible for significant positive GCA effects for KPS under drought-stressed (LM26) and non-stressed (LM48) conditions. Thus, the genetic merit of an individual for KPS can only be maintained under specific soil moisture conditions limiting the usefulness of the above-mentioned parental lines in breeding for increased kernel number under different environmental conditions. Employing several cycles of recurrent selection to increase the allele frequencies of favorable genes for KPS in the current genetic material may be warranted 26 .
The deviation in the expected performance of families as revealed by their SCA's varied greatly among the different sites and water regimes. This confounds the identification and selection of families and individual plants www.nature.com/scientificreports/ for genetic advancement. However, transgressive phenotypes with high SCA effects and involving at least one parent with high GCA effects can be highly heritable 26 . This is important to increase the adaptability of advanced material in water-limited and low-input environments where some parental genotypes perform poorly 59 . Families such as BW141 × LM26 with high SCA effects and at least one parental genotype with high GCA effects for biomass traits and GY may be selected for genetic advancement for drought-prone areas. In addition, reciprocal crosses LM47 × BW152, LM26 × BW140 and LM75 × LM47 with significant reciprocal effects for root biomass and grain yield should also be selected. These genotypes will contribute to soil carbon buildup while concurrently increasing the resilience of wheat productivity in low-input agricultural systems.
Limitations of the study. The plant density per plot was low due to the large number of genotypes evaluated, and the cumbersome process of root sampling through manual excavations further curtailed increasing plot sizes. Data acquired over one year could limit the repeatability of the experiment and the overestimation of variance components and heritability measurements needing multiple environmental analyses. The F 2 generation represents the most segregant populations during the breeding cycle, and individual plants are genetically dissimilar, rendering an increased coefficient of variation of the measured traits. The combining ability analysis was conducted in the F 2 generation representing the maximum genetic variation to discern variance components, and breeding values of parental lines and families, given that hybrid breeding is not the primary objective. The present genetic analysis based on the F 2 generation provided valuable information enabling the selection of promising individuals and families for genetic advancement to develop homozygous and homogeneous lines in the advanced generations.

Conclusion
Significant genetic variation for grain yield, yield-related traits and biomass allocation was observed for the assessed parental genotypes and their families. Additive gene effects conditioned the inheritance of PH, SB, PB and GY under drought, whereas under non-stressed conditions, non-additive gene action was more predominant. Strong maternal and reciprocal genetic effects were recorded for RB across the testing sites under droughtstressed conditions. Genotype LM75 maintained GCA effects in a positive and desirable direction and will be selected for population development in breeding programmes. Early generation selection using PH, SB, PB and GY is recommended to improve drought tolerance by exploiting additive gene action under drought conditions. Higher RB production may be maintained by a positive selection of male and female parents to capture the significant maternal and reciprocal effects detected in this study.

Materials and methods
Plant material, crosses and genetic advancement. Ten genotypes were selected based on their genetic diversity, drought tolerance and ability to produce shoot and root biomass under drought conditions. The pedigree information and the drought sensitivity index (DSI) 15 of the genotypes is summarized in Table 7. Eight of the lines were drought and heat tolerant lines acquired from the CIMMYT drought, and heat nurseries and two lines were local checks with good drought resistance. Crossing blocks of the ten genotypes were established from April to June 2019 at the controlled environment facilities (CEF) at the University of KwaZulu-Natal, South Africa. The parental lines were stagger planted to allow synchronized flowering for emasculation and pollination. A total of 90 families were developed, including 45 direct crosses and 45 reciprocals using a full-diallel mating design. Successful F 1 crosses were harvested from August to October 2019 and bulked to produce F 2 seed in a generation advancement trial from December 2019 to March 2020. Phenotypic evaluation. Field evaluation. One hundred genotypes, including 10 parents and 90 F 2 families, were planted at the Ukulinga Research Farm (29° 40′ S, 30° 24′ E; 806 m above sea level) in a field trial in July 2020. The trial was laid out in a 10 × 10 alpha lattice design with two replications. Fertilizer was broadcasted at a rate of 130 kg N ha −1 , 50 kg P ha −1 and 40 kg K ha −1 before planting to supply plants with adequate nutrition. Test plots were 2 m long and spaced 0.5 m apart, with seven planting stations spaced 30 cm apart. Three seeds www.nature.com/scientificreports/ were planted per planting station and were later thinned to two plants two weeks after germination. Outer field rows were planted with a local cultivar to reduce border effects. An automated drip irrigation system was used to provide water to all plants. The trials were conducted under two water regimes: drought-stressed and nonstressed. The experimental site was covered with a custom-made black plastic mulch to reduce evaporation and stop infiltration of untimely rain water into the soil profile. Drought-stress was imposed by allowing depletion of water in the rooting zone to 35% field capacity at the heading stage of growth in the drought-stressed treatment before rewatering. In the non-stressed treatment, irrigation was continued by maintaining watering at 80% of field capacity until crop maturity. The amount of water in the soil was monitored using a tensiometer (HOBO UX120, Onset, USA) located in each replication to a depth of 60 cm of the rooting zone. The tensiometer readings were used to schedule irrigation in the automated irrigation system for the different water regimes. Standard agronomic practices were kept constant in both water regimes for the duration of the trials according to wheat production guidelines in South Africa 43 . Weather conditions at the site are recorded and reported in 20 .
Greenhouse evaluation. The greenhouse experiment was conducted at the Controlled Environment Facilities at the University of KwaZulu-Natal (29° 37′ S, 30° 24′ E). The experiment was arranged in a 10 × 10 alpha lattice design, with two replications following the pattern of the field experiment. Seven seeds were sown in sterilized composited pine bark growing media in 15L capacity plastic pots with 30 cm and 28 cm height and diameter respectively. Germinated seedlings were thinned to five plants per pot. Fertigation was supplied to the rooting zone using an automated irrigation system to provide sufficient water and fertilizer to the plants. The fertilizer was applied at the same rate as in the field experiment. The experiments were carried out under two water regimes, namely drought-stressed and non-stressed conditions. At the heading stage, stress was imposed in the drought-stressed treatment by allowing the field capacity of the soil to drop to 35%. The non-stressed treatment continued to receive adequate watering to maintain pots at 80% of field capacity until crop maturity. A hand-held moisture probe was used to monitor soil moisture availability in the pots. Insecticides (pyridine azomethine) and fungicides (triazole) were used to control aphids and powdery mildew. Day and night temperatures of 25 °C and 15 °C, respectively and humidity between 45 and 55% were maintained in the greenhouse for the duration of the trial.
Data collection. The following morpho-agronomic traits were measured and recorded on 5 selected plants for every individual plot in the field and pot in the greenhouse: plant height (PH) was measured as the height of the plant from the soil to the tip of the spike using a calibrated meter rule in centimeters (cm); the number of kernels per spike (KPS) was recorded as the number of seeds manually counted from each individual spike. For the two traits, five randomly selected plants were sampled and measured. Shoot biomass (SB) was recorded as the above-ground biomass cut from the base of the plant, including stems and leaves but excluding grain. Root biomass (RB) was recorded as the total root dry matter harvested per genotype per plot. Root samples for each plot were harvested to a depth of 50 cm using a 30 × 30 × 30 cm monolith sampling box. A singe planting station was selected at randomn for harvesting. Large roots were separated manually before washing under running water to remove soil particles. The remaining soil was mixed with water and the suspension was sieved through a 2 mm sieve to collect the fine roots. The fine roots collected from the sieve residue were added and weighed with the large roots. The shoot and root samples were dried at 70 °C for 48 h in an oven drier separately. Dry matter for SB and RB were weighed and expressed in g m −2 . Total plant biomass (PB) was recorded as the total plant dry matter for each genotype in g m −2 . It was calculated by summing up the weight for RB, SB and grain yield (GY), harvested for each genotype. Grain yield was recorded as the total harvested grain per genotype and weighed on a laboratory precision digital scale. The weight of the grain was adjusted to 12.5% moisture content and expressed in g m −2 . Harvest index (HI) expressed in percent was calculated as the ratio of GY to SB, including grain yield as follows: HI = (GY/GY + SB) × 100.

Data analysis. Analysis of variance.
A separate analysis of variance (ANOVA) and Bartlet's test for homogeneity of variance for the two study sites showed significant differences for genotypes and water regimes with homogeneous and comparable variances. Therefore, a combined ANOVA was conducted across the two study sites for the 10 parents and 90 F 2 families using Genstat 18th edition 44 .
Estimation of general and specific combining ability effects. Genetic analysis for a full diallel mating design was computed separately for each test environment using AGD-R statistical software 45 . The GCA and SCA estimates were determined according to 46 Diallel Method I, Model I, following the statistical method below: where: Y ij = phenotypic observation on a cross between the parents i and j, µ = overall mean, g i = GCA effect of parent i, g j = GCA effect of parent j, S ij = SCA effect of a cross between parent i and parent j, r ij = reciprocal effect for the reciprocal crosses between the ith and jth parents, b k = effect of the kth block, e ijk = experimental error due to the environmental effect. The AGD-R software also allowed the further partitioning of the reciprocal effects into maternal (general specific) and non-maternal (specific reciprocal) components.
The relative GCA and SCA ratio was calculated to determine the gene action for each trait using the following formula according to 47 : www.nature.com/scientificreports/ where: σ 2 gca = variance due to GCA and σ 2 sca = variance due to SCA. The broad sense heritability of recorded traits was calculated using the formula below: where H 2 = broad sense heritability, σ 2 g = genetic variance and σ 2 p = phenotypic variance.

Data availability
The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.